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Inverse Spectral Problems for Sturm-Liouville Operators on 
Hedgehog-type Graphs with General Matching Conditions 

V. Yurko 

Abstract. Boundary value problems on hedgehog-type graphs for Sturm-Liouville differen¬ 
tial operators with general matching conditions are studied. We investigate inverse spectral 
problems of recovering the coefficients of the differential equation from the spectral data. For 
this inverse problem we prove a uniqueness theorem and provide a procedure for constructing 
its solution. 

Key words: Hedgehog-type graphs, differential operators, inverse spectral problems 
AMS Classification: 34A55 34B45 34B07 47E05 

1. Introduction. 

1.1. We study an inverse spectral problem for Sturm-Liouville differential operators on so- 
called hedgehog-type graphs with general matching conditions in the interior vertices. Inverse 
spectral problems consist in recovering operators from their spectral characteristics. The main 
results on inverse spectral problems for Sturm-Liouville operators on an interval are presented 
in the monograph [1] and other works. Differential operators on graphs (networks, trees) often 
appear in natural sciences and engineering (see [2] and the references therein). Most of the 
results in this direction are devoted to direct problems of studying properties of the spectrum 
and the root functions for operators on graphs. Inverse spectral problems, because of their 
nonlinearity, are more difficult to investigate, and nowadays there exists only a small number 
of papers in this area. In particular, inverse spectral problems of recovering coefficients of 
differential operators on trees (i.e on graphs without cycles) were solved in [3-4], Inverse 
problems for Sturm-Liouville operators on graphs with a cycle were studied in [5] and other 
papers but only in the case of so-called standard matching conditions. 

In the present paper we consider Sturm-Liouville operators on hedgehog-type graphs with 
generalized matching conditions. This class of matching conditions appears in applications 
and produces new qualitative difficulties in investigating nonlinear inverse coefficient prob¬ 
lems. For studying this class of inverse problems we develop the ideas of the method of 
spectral mappings [6]. We prove a uniqueness theorem for this class of nonlinear inverse 
problems, and provide a constructive procedure for their solution. In order to construct the 
solution, we solve, in particular, an important auxiliary inverse problem for a quasi-periodic 
boundary value problem on the cycle with discontinuity conditions in interior points. The 
obtained results are natural generalizations of the well-known results on inverse problems for 
differential operators on an interval and on graphs with standard matching conditions. 

1.2. Consider a compact graph G in R m with the set of edges £ = {e 0 ,... ,e r }, 
where e 0 is a cycle, £' = (ei,..., e r } are boundary edges. Let {tq,..., v v+ n} be the set of 
vertices, where V = {tq,..., ty}, q 6 et , are boundary vertices, and U = {ry+i,..., u r+ jv} 
are internal vertices, U — £' fl eo . The cycle eo consists of N parts: eo = (J/v=i e r+k , 
e r+ k = [v r +k, xv+AH-i], k — 1, N, v r+ N+i '■= v r+ \. Each boundary edge e 3 , j — l,r has the 
initial point at v 3 , and the end point in the set U. The set £' consists of N groups of edges: 
£' = £fJ.. .U £ n , £fcfleo = v r+ k . Let r k be the number of edges in £ k ; hence r = rq+- ■ ■+r N . 
Denote m 0 = 1, mk — rq + • • • + r k , k = 1,N. Then £ k = {}, j = nrik-i + 1, 777.^, 
Vr+k k = 1 ,N. Thus, the boundary edge e 3 e £k can be viewed as the 

segment e 3 = [v 3 , v r+ k} ■ 

Let Tj be the length of the edge e 3 , j = 1, r + N, and let T := T r+1 + ... + T r+N be 
the length of the cycle e 0 . Put b 0 = 0, b k = T r+1 + ... + T r+k , k — 1, N. Then b N = T. 

Each edge e 3 , j — l,r + N is parameterized by the parameter x 3 6 [0,7}], and x 3 = 0 
corresponds to the vertex v 3 . The whole cycle eo is parameterized by the parameter x 6 
[0, T\, where x — x r+3 + b 3 -\ for x r+3 G [0,T r+J ], j — l,N. 
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An integrable function Y on G may be represented as Y = fa-} - = l r+N , where the 
function y 3 (xj), Xj G [0,7}], is defined on the edge e 3 . The function y(x), x G [0,T] on 
the cycle e 0 is defined by y(x) = y r+ j(x r+3 ), j = 1 ,N. 

Let Q = {qj} J= i ;r+ N be an integrable real-valued function on G ; Q is called the po¬ 
tential. The function q(x), x G [0, T] is defined by q(x) = q r+ j(x r+ j), j = 1 ,N. Denote 
Uj(Y) := yfj{ 0) - hjVji 0), j = 1 ,r + N, U r+N+1 := U r+ 1, where h 3 are real numbers. 
Consider the following differential equation on G : 


-y'-ixj) + q j (x j )y j (xj) = Xy j (x j ), x 3 G [0,7}], j = 1, r + N, (1) 

the functions yj,y'j , j = 1 ,r + N, are absolutely continuous on [0,7}] and satisfy the 
following matching conditions in each internal vertex n M+ i , y = r + 1, r + N : 

^+i(0) = oijyjiTj) for allej G ^_ r+1 , U„ +1 (Y) = ^ /3jZ/'(7}), (2) 

i i 

J y-r -\-1 


yr+N+l • Vr+li ^r+iV+1 • ^r+1, &N+1 ■ 5l i • ^/x—r+1 U e^j, 

where a 3 and b 3 are real numbers, and aj(3j ^ 0. For definiteness, let aj(3j > 0. The 
matching conditions (2) are a generalization of the standard matching conditions (see [5]), 
where a 3 = f3j = 1, h 3 = 0. 

Let us consider the boundary value problem B 0 on G for equation (1) with the matching 
conditions (2) and with the following boundary conditions at the boundary vertices Vi,... ,v r : 

Uj{Y) = 0, j = T~r. 

Denote by A 0 = {A n0 } n >o the eigenvalues (counting with multiplicities) of B 0 . Moreover, 
we also consider the boundary value problems B Vu _ jVp , p — l,r, 1 < v\ < ... v p < r for 
equation (1) with the matching conditions (2) and with the boundary conditions 


Vk{ 0) = 0, k = u 1 ,...,u p , Uj{Y) = 0, j = 1, r, j ± v u „.. , v v . 

Denote by h. Vlj „^ Vp := the eigenvalues (counting with multiplicities) of B Vl ^ Vp 

ft will be shown in Section 2 that an important role for solving inverse problems on graphs 
is played by an auxiliary quasi-periodic boundary value problem on the cycle with discontinu¬ 
ity conditions in interior points. The parameters of this auxiliary problem depend on the pa¬ 
rameters of B 0 . More precisely, let us introduce real numbers 7j, rjj, ( j — 1, N — 1), h,a,/3 : 


a 


r+j 


7j fS ’ ^ Ti^r+j+ij j Ij-Af 1) h h r+ i, 

TV-1 TV—1 TV—1 TV 

OL C^r+TV | 7j | fir+j i ft 11 Tj | fir+j • 

3 = 1 3 =1 3 = 1 j =1 


(3) 


Clearly, a/3 > 0, 7^ >0, j — 1, N — 1. Using these parameters we consider the following 
quasi-periodic discontinuity boundary value problem B on the cycle eo : 


-y"{x) + q(x)y(x) = A y(x), x G (0, T), (4) 

7(0) = ay(T), y\ 0) - hy(0) = /3y'(T), (5) 

y(6j + 0) = 7 jy{bj - 0), y\bj + 0) = 7j"V Oj - 0) + r] 3 y(b 3 - 0), j = 1, N - 1, (6) 

0 < b \ < ... < 6jv —1 &tv = T . 

Let S(x, A) and C(ar, A) be solutions of equation (4) satisfying discontinuity conditions (6) 
and the initial conditions 5(0, A) = C"(0, A) = 0, 5'(0, A) = (7(0, A) = 1. Put <p(x,X) = 
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C(x, A) + hS(x, A). Eigenvalues {A n } n >i of B coincide with zeros of the characteristic 
function 

a(A) = aip{T, A) + (3S'(T, A) - (1 + a/3). (7) 

Put d(X) S(T, A), Q(X) = oap{T, A) — f3S'(T, A). All zeros {z n } n >i of the entire function 
d( A) are simple, i.e. d(z n ) 7 ^ 0, where d(A) := jrd(X). Denote M n = —di(z n )/d(z n ), 
where di(X) := C(T,X). The sequence {M n } n >i is called the Weyl sequence. Put t o n = 

SignQ(^n), = {u n }n> 1 - _ 

We choose and hx one edge G E % from each block £ i jf i = 1, iV, i.e. rrii-i + 1 < £* < 
m,. Denote by £:={£;: k — £ 1 ,..., £at} the set of indices Ci , i = 1, N. Let aj and f3j , 
j — l,r + N, be known a priori. The inverse problem is formulated as follows. 

Inverse problem 1. Given 2 N + r — N spectra A j, j = 0, r , X Ul ,...,u p , p = 2, N, 1 < 
u { < ... < u p < r, Uj G £, and D, construct the potential Q on G and H := [hj\ J= l , r+N ■ 
Obviously, in general it is not possible to recover also all coefficients aq and (dj . Note that 
this inverse problem is a generalization of the classical inverse problems for Sturm-Liouville 
operators on an interval or on graphs. 

Let us formulate the uniqueness theorem for the solution of Inverse Problem 1. For this 
purpose together with q we consider a potential q. Everywhere below if a symbol a denotes 
an object related to q, then a will denote the analogous object related to q. 

Theorem 1 . If Aj = Aj, j = 0 , r, A^...^ = A vi ,...,v p , p = 2,N, 1 < < ... < u p < r, 

Uj G £, and D = D, then Q = Q and H = H. 

We will also provide there a constructive procedure for the solution of Inverse problem 1. 
2. Solution of the inverse problem. 

2.1. Let Sj(xj, A), Cj(xj, A), j = l,r + IV, x 3 G [0,7}], be the solutions of equation (1) 
on the edge e 3 with the initial conditions Sj( 0, A) = C'j(0, A) = 0, A'(0, A) = C}( 0 , A) = 1. 
Put (pj(x, X) — Cj(xj, X) + hjSj(xj, X). Clearly, {(pj(xj, A), Sj(xj, A)) = 1, where (y,z): = 
yz' — y'z is the Wronskian of y and z. 

Let here and below A = p 2 , r := Imp >0, II := {p : r > 0}, Lb := {p : argp G 
[A, 7 r — h]}, 5 G (0, 7 r/ 2 ). For |p| —> oo, p G Lb the following relations hold (see [7]): 

o(A) = t: ]|i . ( 8 ) 

a(A) = 0(e rT ), d{ A) = 0(p~ 1 e rT ), |p| ^ 00 , p G n. (9) 

Fix k = 1, r. Let = {^kj} j= T^+N > be the solution of equation ( 1 ) satisfying ( 2 ) and 
the boundary conditions 

Uj($k) = Sjk, j = l,r, ( 10 ) 

where S ]k is the Kronecker symbol. Denote M k (X) := $fcfc(0,A), k = l,r. The function 
Mfc(A) is called the Weyl function with respect to the boundary vertex v k . Clearly, 

®kk(xk, A) = S k (x k , A) + M k (X)(p k (x k , A), x k e[0,T k ], k = l,r, ( 11 ) 

{tp k (xk, A), $fcfe(x fe , A)) = 1 . ( 12 ) 

Denote M fc L(A) := $ fcj (0,A) , M° kj ( A) := T( ;J (0, A) . Then 

fo, A ) = M kj( X ) S j( x j , A ) + A), Xj G [0, Tj], j = 1, r + IV, k = 1, r. (13) 

In particular, M° fc (A) = 1 , M^ fc (A) = M k (X). Substituting (13) into (2) and (10) we obtain a 

linear algebraic system Lb with respect to M£-( A), z/ = 0,1, j = l,r + lV. The determinant 
A 0 (A) of D k does not depend on k and has the form 

N k 

A 0 (A) = cr(A) fap(A) + II ( 

&=i i</xi<...</x i ( C <Ar i= 1 


( 14 ) 
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^(A) = f[(<Xj<Pj(Tj, A)), ftj( A) = x! ’ a °( A ) = a ( A )’ ai ( A ) = ad ^' ( 15 ) 

j = i a j ( Pjy^ji A) 

We note that the coefficients a 0 (A) and a /il ... /ifc (A) in (14) depend only on S^\Tj, A) and 

Cj U \Tj, A), for j = r + 1, r + IV. We do not need concrete formulae for the other coeffi¬ 
cients a M1 ... Mfc (A). The function A 0 (A) is entire in A of order 1/2 , and its zeros coincide 
with the eigenvalues of the boundary value problem Z? 0 - The function A 0 (A) is called the 
characteristic function for the boundary value problems B 0 . Let A Vli ... );/p (A), p = 1 ,r, 
1 < v\ < ... < u p < r, be the function obtained from Ao(A) by the replacement of 
A) with Sj U \Tj, A) for j = isi,... ,is p , is = 0,1. More precisely, 


N 

A iyi) ... ;i , p (A) = <r vit ...,i/j,(A) ^ao(A) + E E °W —Mfc (A) 

/c=l l</xi<...<Mfc<AT 


k 

x II( E %w + E w)). («o 

i=l ejeS^, j^vi,...,v p ejZEfn, j=vi,—,v P 


<h/l,...,i/p(A) 


n n ( aj Sj(Tj,A)), Si»(A) 

i=l, jj=Vl,...,l/p j = Ul,...,Up 


PjSjiTj, A) 
ajSj(Tj, A) 


(V) 


The function A I/lv .. jI/ (A) is entire in A of order 1/2 , and its zeros coincide with the eigen¬ 
values of the boundary value problem B V1 ^ V . The function A^(A) is called the char¬ 
acteristic function for the boundary value problem B Vli _ >Vp . 

Solving the algebraic system D k we get by Cramer’s rule: M£-(A) = Ajh(A)/A 0 (A) , 
s = 0,1, j = 1 ,r + N, where the determinant A£-(A) is obtained from A 0 (A) by the 
replacement of the column which corresponds to M k -{\) with the column of free terms. In 
particular, 

M fc (A) = -A fc (A)/A 0 (A), k = l~r. (18) 

2.2. It is known (see [8]) that for each hxed j = l,r + N, on the edge ej , there exists 
a fundamental system of solutions of equation (1) {eji(xj, p), 42 ( 4 , p)}, x j £ [0,7}], p G 
II, |p| > p* with the properties: 

1) the functions e^(xj,p), is = 0,1, are continuous for Xj G [0,7}], p G II, \p\ > p* ; 

2) for each Xj G [0,7}], the functions e^(xj,p), is = 0,1, are analytic for Imp > 0, |p| > p* ; 

3) uniformly in x 3 G [0, 7}], the following asymptotical formulae hold 


e ji( x i>p) = (ipY ex v(, i P x j)[ 1 ], e %\ x j,P ) = (—ip) 1 ' exp(—ipxj)[l], p G n, |p| ->• 00 , (19) 


where [1] = 1 + 0(p x ). Fix k = 1 ,r. One has 

$kj(xj, A) = A\. j (p)e j i(x j , p) + A° kj (p)e j 2 (x j , p), x 3 G [0, 7}], j = l,r + N. (20) 

Substituting (20) into (2) and (10) we obtain a linear algebraic system D k with respect to 
A k j(p) , is = 0,1, j = 1 ,r + N. The determinant S(p) of D k does not depend on k, and 
has the form 

, r+N 

S(p) = (h 0 + o(-))p r+Af exp ^ — ip ^ 7}|, (21) 

^ 3 =1 

where ho is the determinant obtained from h(p) by the replacement of e ji(° >P)> e ji( T i>P)i 

€^ 2 ( 0 ,p), Cj^(Tj,p) and hj with 1,0, (—1)", (—1)" and 0, respectively. We assume that 
ho 7 ^ 0. This condition is called the regularity condition for matching. Differential operators on 
G which do not satisfy the regularity condition, possess qualitatively different properties in 
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connection with the formulation and investigation of inverse problems, and are not considered 
in this paper; they require a separate investigation. We note that for classical Kirchhoff’s 
matching conditions we have aq — /3j — 1 , hj = 0 , and the regularity condition is satisfied 
obviously. Solving the algebraic system and using (19)-(21) we get for each fixed x k G 

[ 0 , Tfc) : 

®kk( x k, A) = {ip) v ~ l ex.p(ipx k )[l\, p G n 5 , |p| ->■ oo. ( 22 ) 

In particular, M k ( A) = (ip) _ 1 [l], p G II 5 , |p| —> 00 . Moreover, uniformly in Xk G [0, T*,], 

(p[ u \x k , A) = i((ip)^exp(ipa; fc )[l] + (-ip)" exp(-ipz fc )[l]), p G II, |p| ->■ 00 . (23) 

Using (14), (23), ( 8 ) and (9), by the well-known method, one can obtain the following prop¬ 
erties of the characteristic function Ao(A) and the eigenvalues Ao of the boundary value 
problem Bq . 

1) For pen, \p\ -e 00 : A 0 (A) = 0(exp(r Y^j=i T j ))• 

2 ) There exist h > 0 , C\ > 0 such that |A 0 (A)| > C\ exp(r Tj) for T A h. Hence, 
the eigenvalues A n o = pl Q lie in the domain 0 < r < h. 

3) The number of zeros of A 0 (A) in the rectangle Ag = {p : re [0, h], Rep G 
[£,£ + 1 ]} is bounded with respect to £. 

4) For n —> 00 , p n0 = p ° 0 + O(l/p° 0 ), where A ° 0 = (p ° 0 ) 2 are the eigenvalues of the 
boundary value problem B 0 with Q — 0 and H — 0. 

The characteristic functions A I/lv .. il/p (A) have similar properties. In particular, for p G 
n, |p| ->■ 00 , A I/lv .. )I/p (A) = 0(|p|- p exp(r^+f Tj)). 

Using the properties of the characteristic functions and Hadamard’s factorization theo¬ 
rem, one gets that the specification of the spectrum Ao uniquely determines the characteristic 
function A 0 (A) , i.e. if A 0 = A 0 , then A 0 (A) = A 0 (A). Analogously, if A Vlj _ jVp = A^...^ , 
then A I/ 1 ) ... )l/p (A) = A V 1 ) ... ;i/p (A). The characteristic functions can be constructed as the corre¬ 
sponding infinite products (see [ 1 ] for details). 

2.3. Fix k — l,r , and consider the following auxiliary inverse problem on the edge e k , 
which is called IP(k). 

IP(k). Given two spectra A 0 and A k , construct q k (x k ), x k £ [0, T fc ], and h k ■ 

In IP(k) we construct the potential only on the edge e k , but the spectra bring a global 
information from the whole graph. In other words, IP(k) is not a local inverse problem related 
to the edge e k . Let us prove the uniqueness theorem for the solution of IP(k). 

Theorem 2 . Fix k = l,r. If A 0 = A 0 and A k = A k , then q k (x k ) = q k (x k ) , a.e. on 
[0, Tfi, and h k — h k . Thus, the specification of two spectra A 0 and A k uniquely determines 
the potential q k on the edge e k , and the coefficient h k . 

Proof. Since A 0 = A 0 , A* = A k , it follows that A 0 (A) = A 0 (A), A*,(A) = Afc(A), and 
according to (18), 

M k ( A) = M k { A). (24) 

Consider the functions 

Pi s (x k ,X) = ( _ l ) s_1 (ffk{x k , X)^ ( k f s \x k , A) - fi k ~ s \x k , \)<& kk (x k , A) j, s = 1 , 2 . (25) 

Using (12) we calculate 


Vk(Xk, A) = Pn(Xk, A )ih{Xk, A) + P, 2 (x k , X)ip' k (x k , A). 


(26) 


It follows from (23), (24) and (26) that 


Pi s ( x ky A) = A ls + 0(p '), peH f , |p| -> 00 , x t e (0, T k ■ 


(27) 
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According to (11) and (25), 

P?s(x k ,\) = (-lr 1 ((W, X)st S \x k , A) - <ft S \x k , X)S k (x k , A)) 

+(M k ( A) - Mfc(A))^ fc (x fc , A)<p|j. 2 _s) (a; fc , A)). 

It follows from (24) that for each fixed x k , the functions Pf s (x k , A) are entire in A of 
order 1/2. Together with (27) this yields Pft (x k , A) = 1, Pf 2 (x k , A) = 0. Substituting 
these relations into (26) we get ip k (x k ,A) = (p k (x k ,A) for all 27 , and A, and consequently, 
q k {x k ) — q k (x k ) a.e. on [0,7/], h k = h k . Theorem 2 is proved. 

Using the method of spectral mappings [ 6 ] for the Sturm-Liouville operator on the edge 
e k one can get a constructive procedure for finding q k and h k . 

Now we study the following auxiliary inverse problem on the cycle e 0 , which is called 
IP(0). Consider the boundary value problem B of the form (4)-( 6 ), where the parameters 
of B 0 are defined by (3), and a, (3 are known. 

IP(0). Given a(A),7(A) and fl, construct q(x), x G [0, T], h, 7 ^ and r}j,j — l,N — l. 
This inverse problem is a generalization of the classical periodic inverse problem. Moreover, 
for the standard matching conditions ( otj — (3j — 1, hj — 0 ), IP(0) coincides with the classical 
periodic inverse problem. 

This inverse problem IP(0) was solved in [7], where the following theorem is established. 

Theorem 3. The specification a(A),7(A) and U uniquely determines q(x),h,'fj and 
rjj , j — 1, N — 1. The solution of IP(0) can be fowid by the following algorithm. 

Algorithm 1. 

1) Construct D( A) = a(A) + (1 + a(3). 

2) Find zeros {z n } n > i of the entire function 7(A). 

3) Calculate Q(z n ) via Q(z n ) = u n ^JD 2 (z n ) — 4afi. 

4) Construct di(z n ) by dfiz n ) = (D(z n ) + Q(z n ))/(2a). 

5) Find d(z n ). 

6 ) Calculate the Weyl sequence {M„,} n >i via M n = —di(z n )/d(z n ). 

7) From the given data {z n ,M n } n >i construct q(x),'yj,r}j , j = 1, A — 1, by solving the 
inverse Dirichlet problem with discontinuities inside the interval (see [9]). 

8 ) Find S(T, A), S'(T, A) and C(T, A). 

9) Calculate h, using (7). 

2.4. Let us go on to the solution of Inverse problem 1. Firstly, we give the proof of 
Theorem 1 . Assume that A k = A*, k — 0, r, A„ lt _ jVp = A„ lt _ jVp , p = 2 , N, 1 < iq < ... < 
v p < r, iq G /, and = U. Then one has 

A fc (A) = A fc (A), k = Opr, 

A vl ^ p (X) = A J/lj ... ) 1 /p (A), p — 2,N, l <v 1 < ... <v p <r, v 3 e £. 

Moreover, according to (3), 7 j = 7 j , j = 1, A — 1, and a = a, ft — ft. Using Theorem 2, 
we get q k (x k ) = qk(x k ) a.e. on [0,7/] and h k = h k , k = 1 , r, and consequently, 

Cfc( 27 , A) = C fc (a; fc , A), 5/(a; fc , A) = Pjt(x fe , A), <p fc (a; fc , A) = ^ fc (x fc , A), k = 1, r. (28) 

By virtue of (15), (17) and (28) one has 

<r(A) = flf(A), .^(A) = .^(A), U,(A)^U,(A), n°(A) = «J(A), j = iff. 

Using (14) and (16), we obtain, in particular, a 0 (A) = a(A) , ai(A) = ai(A). In view of (15), 
this yields a(A) = a(A), 7(A) = 7(A). It follows from Theorem 3 that q k (x k ) = q k (x k ) a.e. 
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on [0, Tk], k — r + 1, r + N, and h = h, rjj = fjj, j — 1, N — 1. Taking (3) into account, we 
get 77 = 77. Theorem 1 is proved. 

The solution of Inverse problem 1 can be constructed by the following algorithm. 
Algorithm 2. Given A*,, k = Opr, A V1 ^ V , p — 2, N, 1 < iq < ... < u p < r, Uj e £ t 
and O. 

1) Construct A k (A) and A^.^A). 

2) Calculate 7 j, j = 1, N — 1, a and /?, using (3). 

3) For each fixed k = 1 ,r, solve the inverse problem IP(k) and find q k (xk), %k £ [0, T k ] on 
the edge e k and h k . 

4) For each fixed k — l,r, construct C k (x k , A), S k (x k , A) and (p k (x k , A), e [0, T fc ], 

5) Calculate a(A) and d(A), using (14), (15) and (16). 

6) From the given a(A),d(A) and 11, construct q k (x k ), [0,7*.], k = r + l,r + N, h and 

Vj , j = 1,77-1- 

7) Find 77, using (3). 
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